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Total-angular-momentum (TAM) waves provide a set of basis functions for scalar, vector, and 
tensor fields that can be used in place of plane waves and that reflect the rotational symmetry of 
the spherical sky. Here we discuss three-point correlation functions, or bispectra in harmonic space, 
for scalar, vector, and tensor fields in terms of TAM waves. The Wigner-Eckart theorem dictates 
that the expectation value, assuming statistical isotropy, of the product of three TAM waves is 
the product of a Clebsch-Gordan coefficient (or Wigner-3j symbol) times a function only of the 
Qf) , total-angular-momentum quantum numbers. Here we show how this works, and we provide explicit 

expressions relating the bispectra for TAM waves in terms of the more commonly used Fourier-space 
bispectra. This formalism will be useful to simplify calculations of projections of three-dimensional 
fSJ , bispectra onto the spherical sky. 

in 



O 



Oh 



C3 



> 

O 



- 1—1 

x 



I. INTRODUCTION 



Cosmological measurements have over the past several decades made the notion of a period of inflationary expansion 
in the early Universe particularly appealing. The principal aim of early-Universe cosmology has thus become the 
elucidation of the new physics responsible for inflation. While the simplest single-field slow-roll (SFSR) inflation 
' models predict primordial perturbations to be very nearly Gaussian [l[ , they do require some nonvanishing departures 
from Gaussianity @- Moreover, just about any extension of SFSR models, or embedding of toy SFSR models into 
more realistic models or alternatives or additions to inflation Q lead to larger departures from Gaussianity. 

i - "! There has thus been ever-growing attention focused on the search for non-Gaussianity. 

The vast majority of the literature on departures from Gaussianity focuses on the bispectrum (the three-point 
correlation function in Fourier space) for the primordial curvature, a scalar quantity Q- However, inflation also 
predicts primordial gravitational waves (tensor metric perturbations) Q, and some inflationary models involve the 
introduction of vector fields [1] . There are also discussions of primordial magnetic fields ■ Just as there may arise 
non-Gaussianity in the primordial scalar perturbation, there may also be non-Gaussianity in these vector fields, in 
magnetic fields [HI, [3, and in the gravitational- wave background [T^ - fToT ]. Additional work on non-Gaussianities with 
vector and scalar fields can be found in Refs. [l6j . 

Non-Gaussianity in the primordial scalar perturbation is most commonly parametrized in terms of a bispectrum, 
the expectation value for the product of three Fourier modes of wavevectors kx, k2, and k3, of the fields under 
\Q , consideration. If the fields under consideration include vector or tensor fields, then the bispectrum will also depend 
on some contractions of the polarization vectors/tensors for these fields. 

Many cosmological observations, however, are performed on a spherical sky, and many obscrvables (e.g., CMB tem- 
£N) ■ perature/polarization) depend only on an angular position on the sky. While others may depend on three-dimensional 
positions, astronomical measurements discriminate between the two angular coordinates and a radial coordinate. 
Comparison of measured quantities with theoretical Fourier-space bispectra therefore necessarily require projection of 
the three-dimensional Fourier-space bispectra onto the two-dimensional spherical sky. The bispectra projected onto 
the spherical sky depend on the three multipoles l\, I2, and Z3, associated with the three spherical-harmonic coeffi- 
cients being correlated. Such projections necessarily involve the complications associated with spherical harmonics, 
Clebsch-Gordan coefficients, and Wigncr-3j symbols, and sometimes even Wigner-6j and Wigner-9j symbols. 

In a recent paper (l7| , we developed a total-angular-momentum (TAM) formalism to provide sets of basis functions 
for scalar, vector, and tensor fields in three spatial dimensions that reflect the rotational symmetry of the Universe 
about any given point (taken to be our location). The purpose of these TAM waves is to incorporate the rotational 
symmetry of the sky from the start, rather than start with plane waves which are later projected onto the sky. These 
TAM waves generalize the more familiar tensor spherical harmonics fl8l - [20| , which provide a basis for tensor functions 



on the two-sphere, to a basis for functions in three-dimensional Euclidean space. The TAM waves of Ref. [17[ are 
also provided for scalar and vector fields, and for all five components of a traceless tensor, not just the two transverse 
components. 

Each TAM wave is labeled by a wavenumber k, total-angular-momentum quantum numbers J, and azimuthal 
quantum number M . There are three sets of TAM waves for each kJM for vector fields (to reflect the three components 
of a vector field) and five sets of TAM waves for each kJM for a symmetric traceless rank-2 tensor. We provided in 
Ref. [I?]] several different sets of bases, for a given JM, for these three vector TAM waves and five tensor TAM waves. 
In the orbital-angular-momentum (OAM) basis, the three vector basis functions for a given JM have orbital angular 
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momentum I = J — 1, J, J + l, and the five tensor basis functions have I = J — 2, J — 1, J, J+l,J + 2. In the second 
basis, the three vector basis functions for a given JM include a longitudinal (L) mode and two transverse modes (E 
and B) of opposite parity. The five tensor functions in this basis include a longitudinal (L) mode, two vector modes 
(VE and VB) of opposite parity, and two transverse-tensor modes (TE and TB) of opposite parity. The third basis 
represents vector fields of a given JM in terms of states of helicity A = 0, ±1 and tensor fields in terms of modes of 
helicity A = 0, ±1,±2. The most general scalar, vector, or tensor field can then be written in terms of these TAM 
waves, rather than Fourier modes. Spherical-sky observables are then obtained more naturally from these TAM waves 
than from Fourier waves. 

In this paper we calculate the bispectra for TAM waves in terms of the more commonly used Fourier-space bispectra. 
As we will see, the Wigner-Eckart theorem guarantees that the expectation value of the product of three TAM-wave 
coefficients depends on M 1; M 2 , and M 3 only through the Clebsch-Gordan coefficient (or equivalently Wigner-3j 
symbol or Gaunt integral); there may then be some dependence on J\, J 2 , and J3 in a prefactor in addition to that 
in the Clebsch-Gordan coefficient. The specific form for the prefactor will depend on whether the bispectrum is 
for scalar, vector, and/or tensor fields as well as on the tensorial nature of the bispectrum (i.e., how the indices on 
the polarization vector/tensors are contracted). The principal results of this paper are thus the specific forms for 
these prefactors for a variety of bispectra. These results will facilitate the calculation of observables, such as angular 
bispectra, for theories that involve such correlations, particularly those that involve vector and/or tensor fields. 

We begin in Section [TT] with a brief review of the total-angular-momentum wave (TAM) formalism, and we use 
the notation and conventions of Ref. [l?} throughout. Section IPTE! presents the bispectra for three scalar TAM waves. 
Section HVl follows to calculate bispectra involving two transverse- vector TAM waves, in the E/B basis, and one 
scalar TAM wave. Section [V] presents bispectra of two transverse-traceless tensor TAM waves plus one scalar TAM 
wave. Section IVT1 then deals with the correlation of one symmetric traceless tensor TAM wave with two scalar TAM 
waves, where the traceless tensor can be a transverse-traceless, vectorial, or longitudinal. The same results apply 
also to correlation of one traceless tensor TAM wave and two longitudinal vector TAM waves. Section IVlII provides 
concluding remarks and provides a list of the equation numbers for the central results that relate the Fourier-space 
and TAM-wave three-point functions. Appendix [^presents helicity-basis overlap integrals that are used earlier in the 
paper. We discuss a semi-classical picture to understand the structure of the TAM-wave bispectra in Appendix IB! 



II. REVIEW OF TOTAL ANGULAR MOMENTUM WAVES 



Here we very briefly re-introduce total-angular-momentum (TAM) waves for scalar, vector, and tensor fields. We 
follow throughout the convention and notation of Ref. |17| and refer the reader there for more details. 

TAM waves are eigenfunctions of eigenvalue —k 2 of the Laplacian that are also eigenstates of quantum numbers 
J and M of total angular momentum squared and its third component. The scalar TAM waves are denoted by 



*5 JAr) (x). TAM vector waves of orbital angular momentum I are specified by ^/ /jU N (x), where a = {x,y,z} is the 



vector index and I = J — 1,J, J+l. TAM tensor waves of orbital angular momentum I are specified by \E'|j M j ab (x), 
where a and b are the tensor indices and I = J — 2, J — 1,J, J+l,J + 2. The scalar, vector, and tensor TAM waves 
are distinguished by the number of indices. 

These TAM waves satisfy the completeness relations, 
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These relations imply that an arbitrary scalar, vector, or tensor field can be expanded, respectively, 
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in terms of TAM waves and expansion coefficients, 
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By using the plane- wave expansions, 
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and orthogonality of the TAM waves, 
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we can also write the expansion coefficients for scalar, vector and tensor waves as 
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Similarly, if we choose to decompose the vector field into L/E/B waves and the tensor fields into L/VE/VB/TE/TB 
waves, which we label by a, then we also have 
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(JM) 
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Finally, we can write the L/E/B vector TAM waves in terms of the L/E/B vector spherical harmonics as 



*UM)a( x ) = Jj(kr)Y ( % )a (n) 



(JM) 
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and the L/VE/VB/TE/TB tensor TAM waves in terms of the L/VE/VB/TE/TB tensor spherical harmonics as 
*&M)a 6 ( x ) =-\ (Mkr) + 3j'j(kr)) Y { L JM)ab (h) - ^3J(J + l)fj(kr)Y™ )ab (n) 



1 3(J + 2)\jj(kr) vTE 



(J -2)! (kr) 2 (JM '> 



*Ko«*(*) = ^J(J+l)fj{kr)Y ( L JM)ab (h) (3.;(kr) + 2j'j(kr) + 2fj(kr)) Y^, E M)ab (h) 
V(J-l)(J + 2) (jj(kr) + 2 J ^Pj Y { T jf I)ab (&), 

- \ {-jAkr)+mr)+^Akr)+Q 3 j^j Y? E M)ab {n), 

*Jj b mU*) = - i (fjM - ^r 1 ) r ( ^ )a6 (fi) - iv / FWT2)^^,.(fi) ! 

= - 'vTW^F ( S )aJ (n) - i (j'Akr) + 2^M) F ( ^ )a6 (fi), (22) 

where fj(x) = (j j(x)/x)' . For later use, we define additional radial functions for vector TAM waves with a,/3 = E, L 
by 

= -iE^w^W' (23) 

and for tensor TAM waves with a, j3 = L, VE, TE by 

*(&)o6M 55 - EiS^^^W. ( 24 ) 

and for tensor TAM waves with a, (3 = TB, VB by 

^M )ab « = "iE^^^W- ( 25 ) 

The precise forms for jj*^ can be read off Eq. (|2"Tj) and for jjf from Eq. (|2"2"j) . Some of these radial profiles have 
appeared in previous full-sky studies in cosmology. For example in the line-of-sight approach to CMB polarization, the 
full-sky EE and BB power spectra, Eq.(29) of Ref. [I(J), contain radial kernels identical to the functions jj t E ' T (%) 
and jjt B ' TB \%) defined here. 

A. From TAM-bispectrum to angular bispectrum 

Throughout this paper, we will calculate TAM bispectra of the form, 

.VliWi J9M9 JiMz / 1 ^ ' 



Ji Mi 1 J 2 M 2 ^J 3 M 3 j i 

where X, Y , and Z represent generic scalar, vector and tensor fields, and <Xi represents the particular mode of interest 
(L/E/B for vector fields and L/VE/VB/TE/TB for tensor fields). The TAM-wave bispectra above are related to 
the angular bispectra for various observables as follows. Suppose we observe some projection of a three-dimensional 
field -X"(x). The spherical-harmonic coefficients for the projection can be written in terms of the Fourier components 
X(k) of the field as 

/d 3 k f k 2 dk 

^ 5x (fc)X(k)3C*(k) = J j^9x(k)X^ p (27) 
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in terms of a transfer function gx(k) which, assuming statistical isotropy, is a function only of the wavevector 
magnitude k. Here we suppress tensor indices, as the end result (as will be spelled out more clearly in the rest of 
the paper) is the same for scalar, vector and tensor fields, and yf m gcnerically represents the corresponding a modes 
of scalar, vector, and tensor spherical harmonics. The observable angular bispectrum is then related to the TAM 
bispectrum as 

/ ai,X a 2 .Y a 3 .Z \ [ k x dk\ f k%dk 2 n f k%dk 3 \/v«i,fav«2,fa7a3,te\ (oe\ 

("jiM^'mMm,/ = J J^9x(ki) J j^9Y(k2) J 1^9z{kz) {X^fc Zj 3 ' Ms j . (28) 

For example, if we take X(x) to be Bardeen's curvature perturbation $(x), and the af m to be CMB-temperature 
spherical-harmonic coefficients, then the appropriate transfer function is 

g T (k)=4Wg Tl (k), (29) 

in terms of the radiation transfer function gri(k) [2lJ. If the af m are taken to be the density of galaxies on the sky, 
then the transfer function is 

g g (k) = Wb g (^^) / dzW g (z)D(z)Mkd A (z)}, (30) 

with galaxy bias b g , matter transfer function T(k)„ linear growth factor D(z), comoving angular=diameter distance 
cIa(z), and redshift distribution W g {z) of galaxies normalized to / dzW g (z) = 1. Here, Hq and fl m are, respectively, 
the Hubble and matter-density parameters at present. 



III. SCALAR BISPECTRUM 



Before turning to vector and tensor fields, for which the TAM formalism provides the most substantial advantage, 
we begin by way of introduction with scalar bispectra in the TAM formalism. 

The bispectrum B sss (ki, k2, k 3 ) for the scalar field is usually defined in terms of Fourier modes as 

(0(ki)0(k 2 )0(k 3 )) = (2^) 3 Mkx + k 2 + k 3 )B sss [k!, k 2 , k 3 ), (31) 

where the Dirac delta function arises as a consequence of statistical homogeneity, and the dependence of the bispectrum 
only on the magnitudes ki, fe, and k 3 as a consequence of statistical isotropy. The subscript 'sss' is used to denote 
the bispectrum for three scalar fields, to distinguish it from the bispectra to be discussed below that involve vector 
and/or tensor fields. 

If the scalar bispectrum is defined as in Eq. (|3T|) . then the expectation value for the product of three TAM- wave 
coefficients are given by 

{^UiAXm 2) ^m 3) ) = J d 2 U J d 2 k 2 J d 2 k 3 (0(k 1 )^k 2 )^k3))r ( : /iA/i) (kor ( * j2 M 2 )(k 2 )r ( * j3M3) (k 3 ). (32) 

We then substitute Eq. ([31]) and expand the Dirac delta function, using the plane- wave expansion, Eq. (fTO)) . as 
(2^) 3 fo(ki+k 2 +k 3 ) = y*d 3 a;e- ?;klX e- lk2X e- lk3X 

= / ^S(47r) 3 (-0 ,1+ ' a+,8 ii I (fcir)ii 3 (fcar)i Js (A ! 3r) 

l-irrii 

to eliminate the k^ integrals. We can then write the TAM- wave bispectrum as 

(^mAU.Aw) = (^f(-i) Jl+j2+j3 B sss (ki,k 2 ,k 3 ) J d 3 x*f}; Mi) (x)^5; M2) (x)^5; M3) (x), (34) 

in terms of an overlap of three TAM waves. Below we will see that the TAM-wave three-point functions will always 
depend, even with vector and tensor fields, on overlaps of three TAM waves. Since the scalar TAM waves are 
*jm(x) = jj(kr)Y {JM) (x), the overlap is simply, 

/ ^^M.W^ifcW^t^W^Sfs^Waft-^.y, (35) 
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where 

W a l, f ,2^ ({iWi ™ m / (2? 1 + 1)(2Z 2 + 1)(2Z 3 + 1) / h h h\(hhh\ 
Gm^rnz -J d nY (ilTOl) (n)% m2 ) (n)F ((3ro3) (n) - y [rm m 2 m 3 ){o ) ' 

(36) 

is the Gaunt integral, and 

/•OO 

Ji 1 i 2 i 3 (ki,k 2 ,k 3 ) = / r 2 drji 1 (k 1 r)ji 2 (k2r)ji 3 (k 3 r). (37) 
Jo 

Therefore, 

(^MO^M^^Ms)) = ( 4 ^) 3 (-0 J1 + J2 + J3 ^a} 1 m'm 3 -Wafo, fc 2, k 3 )B sss (k u k 2 , fcs). (38) 

We have thus evaluated the TAM-wave three-point correlation function in terms of the overlap of three TAM waves. 
We then arrive at our final result, Eq. (|38[) . which provides a simple expression for the TAM-wave three-point function 
in terms of the Fourier-space bispectrum, by evaluating that overlap. The factors connecting the Fourier-space and 
TAM-wave three-point functions are the easily evaluated and familiar Gaunt integral, and the radial integral, Eq. (|37l) . 
The Gaunt integral is proportional to a Wigner-3j symbol and the Clebsch-Gordan coefficient. Its presence here is 
a result of the Wigner-Eckart theorem, and we will see in the Sections below that the expectation value for the 
product of any three TAM waves — not just scalar TAM waves — will be proportional to a Clebsch-Gordan coefficient 
times a function only of J±, J 2 , and J 3 . We choose to write the results in terms of the Gaunt factor, rather than 
Clebsch-Gordan coefficients or Wigner-3j symbols simply because it provides for more com pac t expressions. 

The radial integral in Eq. ([3"5)l can be evaluated numerically using standard techniques [22H24| . but often an end 
result, which requires an additional integral over the ki, can be obtained more quickly by changing the orders of those 
ki integrals and the radial integral in Eq. Q37p. 

Below we will generalize Eq. ([58)1 to three-point functions involving vector and tensor fields as well. 



IV. TWO TRANSVERSE VECTORS AND A SCALAR 

We now proceed with the analogous calculation for the three-point correlation function that involves two transverse- 
vector fields and a scalar. The transverse vector field V a (x) satisfies V a V a — 0, or in Fourier space, k a V a (k) = 0. 
Such a vector- vector-scalar correlation arises, for example, in models in which a magnetic field (the transverse- vector 
field) is produced during inflation and thus correlated with the inflaton field (the scalar field) [n| . 

Our goal is to calculate the TAM-wave bispectra, 



for {a\,a 2 } = {E,B} that correspond to a given Fourier-space bispectrum. 

A. Fourier-space bispectra 

We begin by discussing the Fourier-space bispectrum. Each Fourier component of the transverse-vector field can 
be written 

V a (k)= £e a A (k)V> A (k), (40) 

in terms of two helicity-basis polarization vectors i\(k) and hclicity Fourier amplitudes V\(k). The helicity basis 
vectors can be written as e±(k) = £i(k) ± ik x £i(k) /V2, where i\ is any unit vector in the plane perpendicular 
to k. The most general Fourier-space vector-vector-scalar three-point function can therefore be written, 



Ai A2 



^(k 1 )F 6 (k 2 )^(k 3 )) = > el 1 (k 1 ) £ l 2 (k 2 )(^ 1 (k 1 )Vo, 2 (k 2 )^k 3 )). (41) 
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Given the orthonormality of the polarization vectors, this relation can then be inverted to provide the three-point 
correlation, 

(VA I (k 1 )V Aa (k a )0(k 3 )) = £ X(ki)4 2 (k 2 ) (K(k!)n(k 2 )0(k 3 )) , (42) 

for the helicity amplitudes. 

Now consider the tensor structure of the three-point function on the right side of Eq. (|42|) . By momentum conser- 
vation, it most generally depends on and k%- Among all structures one can possibly construct, g ab , k^k\, k^k\, 
k^k\ and k 2 k\, terms proportional to k\ or k 2 have no contribution to Eq. (|42j) under helicity basis, and hence to the 
bispectrum through Eq. ([SJ). The most general parity-conserving three-point function can then be written, 

V a {ln)V b {k 2 )^ 3 )) = (2 7 r) 3 fo(k 1 + k 2 + kg) \g ab B^ s {k u k 2 , k 3 ) + k 2 ,k 3 )] , (43) 



in terms of two vector-vector-scalar bispectra BTOs(fci,fe]fe) and Bvv a (ki,k2,ka). Two different bispectra arise 
because there are two components for each polarization vector. 

To demonstrate our formalism, we take the example of non-zero B^ a (k±, k2, ks), which may arise from a local 
three-point interaction of the form </>(x)|V(x)| 2 . The other bispectrum B$ s {k\, k%, kz), which can arise as the vector- 
vector-scalar correlation between the vector eletromagnetic potential and the inflaton from a (f>F IJiV F l '' v coupling 
is left to future work. 

Therefore, the Fourier-space bispectrum we will consider here will be of the form, 

y a (k 1 )V> b (k 2 )0(k 3 )\ = (2^) 3 fe(k 1 +k 2 + k 3 )g ab B vvs (k 1 ,k 2l k 3 ). (44) 



B. TAM-wave bispectra 

We now calculate the three-point correlation function, 



tftft)*8£>*&*,)) = / d 2 k ld ^k 3 ^M l)a (k!)y ( -4 )6 (k 2 )F ( ; 3M3) (k 3 ) (^(k 1 )^(k 2 )^k 3 )) , ( 45 ) 

for the E/B TAM-wave coefficients, where ai, a.% = {E, B}, and we have used Eq. (fT9|). The first step is to expand 
the Dirac delta function in Eq. (|4"4"|) . and using the identity, 



J'M' 

we find that 



E (-iy'M^Yfr^in) / d^Yft^&Y^^ik) = (-i) J fe f)o (x)l * , (46) 



(^ft)^Sft)^v.)) = (^?H) Jl+M B vvs (kuk 2 ,k 3 ) J d 3 ^-^-(x)^^ a (x)^ 3 * M3) (x). (47) 

That is, the bispectrum of three TAM coefficients is proportional to the overlap of three TAM waves. 

These overlap integrals can be evaluated with the projection of vector TAM waves in Eq. (|2"Tj) . The angular parts of 
the overlap integral can be evaluated by writing the L/E/B vector spherical harmonics in terms of helicity spherical 
harmonics as 

(JM)a = (^JM)o ± iY (JM)a) A^! Y (!lM)a = Y (JM)a- (48) 



1 If parity violation is allowed, there may be additional terms constructed with the antisymmetric tensor e a j, c , but we restrict our attention 
to parity-preserving terms. 
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We can then use the helicity-harmonic integrals in Appendix I A II to obtain the spherical-harmonic overlap integrals, 



d2fi ^; a M l) ( fi ) y ( L / 2 A/ 2 ) Q ( fi ) r (^M3)(fi) 

rf 2 ^5M 1 )( ft ) y (5 2 M 2 )a(fi)^(^M3)(fi) 



t Ji J 2 J3 



&M l M 2 M 2 > 

d2ft ^M 1 )( fi ) r (iM 2 )a(fi)^(J3M 3 )(fi) 



l+(-l) 



Jl+J2+J3 C J M (JiU2,-l|J 3 0) 

2 y Ml M 2 M 3 ( Jl0 j 2 o|j 3 o) ' 



* M1M2M3 ' 



( jq-h, -i\j 3 o) 

<JiOJ 2 0|J 3 0) ; 



; \j\m, ) ( fi ) y (./ 2 A/ 2 )a ( fl ) y (Wi 3 ) (") 

and we also have that 

F ( ^ iMi)a (n)y ( ^(n)=0. 
From these we obtain that the nonzero TAM-wavc overlap integrals for Ji + J 2 + J3 = even arc 



I 



^^Sm i)0 ( x )*& )( x )*& 3 m 3 )W 



= g 



> Ml M 2 M 3 (J l0 J 2 0|J 3 0) 

(JilJ 2 ,-l|J 3 0) 

(J!0J 2 0|J30) 



-JjlJ 2 J 3 { k l-, k 2, k 3 ), 



(49) 



(50) 



(51) 



J1J2J3 

Ah Ah M3 



r * dr 3j'v ( k i r )j)h!v' ( k 2r)jj 3 (for) 



■ (E,E) 



^drjjuv' ( k i r )3jt!v' (for)jj 3 (for) 



(EX) 



(52) 



and when J\ + J 2 + J3 = odd, 

/ ^A a to*!**) to 



||^ ga 3 / - 2 drj Jx for) i£f J (k 2 r) jj3tV (k 3 r) 



(53) 



if Ji + J2 + J3 = odd, as the Wigner-3j symbol 



vanishes. In every case where a Gaunt integral appears 



Here, wc have used the radial functions defined in Eq. (|23|) . Note that strictly speaking, the Gaunt integral vanishes 

Ji J2 J3 


with Ji + J 2 + J 3 = odd, though, it appears divided by a Clcbsch-Gordan coefficient (Ji0J 2 0| J3O) which also contains 
the same Wigner-3j symbol, and thus cancels that in the numerator. We choose to write things in this slightly 
unconventional way to enable compact, but still unambiguous, expressions. 

Eq. (|47|) . together with the overlap integrals in Eqs. (|5ip . (f52|) . and (|53|) . provide the £?-i?-scalar, £'-£ , -scalar, and 
.B-E'-scalar TAM-wave bispectra in terms of B vvs (ki, k%, fo), the Fourier space vector-vector-scalar bispectrum. 



V. TWO TRANSVERSE-TRACELESS TENSORS AND ONE SCALAR 

We now consider a three-point correlation involving two transverse-traceless tensors and one scalar. Such a corre- 
lation arises, for example, in inflation if the tensor modes are gravitational waves and the scalar field is the curvature 
perturbation [l3j . 

Our goal here will be to calculate the TAM-wave bispectra, 

frpaiM rpa 2 ,k 2 ,k 3 \ (ZA\ 

Y (JiAfi) J {J 2 M 2 ) ( P(J 3 M 3 ) J ' 

for {ai,a 2 } = {TE,TB} that correspond to a given Fourier-space bispectrum. Here we will take the Fourier-space 
tensor-tensor-scalar three-point function to be of the form 

(f «*(ki)f cd (k 2 )0(k 3 )) = (27T) 3 6 D (k 1 + k 2 + k3) 3 ac 5 M -Btt,(fci, fo, fo), (55) 

that arises in single-field slow-roll inflation in terms of a tensor-tensor-scalar bispectrum B Us {ki, fc 2 , k 3 ). The 
tensor fields T ah are transverse, V a T a& (x) = 0, and traceless, g ab T ab (x) = 0. 
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We calculate the three-point correlation function of TE/TB TAM coefficients directly from Eq. ([20)) : 

( T (7;m\) T (Xm 2 )^3M 3 )) = / rf 2 kAd 2 k 3 F-^ (56) 
where ax,o>2 = {L,E, B}. Expanding the Dirac delta function in Eq. (|55p with Eq. (|33p and using the identity, 

E (-*) J '^(fcO^M')( fi ) / rf 2 k^M )a6 (k)^(j'M')(k) = H) J [^j k M)ab (x)} * , (57) 

7/ n r/ « 



J'M 

we find that 



lAw) = (4 7 r) 3 (- l )^ + ^ + ^^(fc 1 ,fc 2 ,fc 3 ) / d 3 x*- ; ^-^(x)*^^ ab (x)*^ M3) (x). (58) 



1 (JiM 1 ) 1 (J 2 M 2 ) 

That is, the bispectrum of three TAM-wave coefficients is again proportional to the overlap of three TAM waves. 
Our task now is to evaluate the overlap integrals, 

l3„ir,ai,ki I \ T a 2 ,k 2 ,ab 



^TjIm^Wj,m 2 ] °(*)*(,/3M3)( x )' (59) 

for a,\,a.i = TE,TB. With the the decomposition in Eq. (94) of Ref. [17j , we only have to calculate the overlap of 
three tensor spherical harmonics 

"d 2 nJ^ cd y- Mi)o6 (n)y- M2)cd (fi)F (j3 M3)(fi), (60) 

for ot\, a.2 = L, VE/VB, TE/TB. Here the constant tensor I ab ' cd is the identity tensor acting on symmetric traceless 
tensors. In Appendix I A 21 we calculate these overlaps in the helicity basis. We then use the relations, 



Y ±2 _ 1 ( Y TE +iY TB \ Y^ - — (y ve + iY VB \ Y a -Y L (M) 
I (JM)ab ~ ^ V ( JM ) ab ll (JM)abJ ' 1 (JM)ab ~ V ( JM ) ab L1 (JM)ab) > 1 (JM)ab ~ 1 (JM)ab' 

to find that 

' d 2 flY^{m { \M 2 )ab{W{JsM 3 ){Z) = gtiXM 3 , (62) 
^ v VBab (A\vVB - / r fi^vVE,ab ,~^TrVE 



3 nF ( ^(fi)F ( ^ M2)ai) (fi)F( j3 M3)(fi) = j ^nF ( ^X)( fi )^ 2)a 6(fi)>(73M 3 )(fi) 

2 y Ml M 2 M 3 ( j l0 j 2 0|J 3 0) ' lbdj 

J2- yVB,ab V VE (~\ Vi _ 1 - (- 1 ) Jl+ ' /2 + ' 73 r ,J 1 J 2 J3 ( J 1 1J 2,-1|J30) , R .x 
d ny (J 1 A/ 1 ) y (J2Af 2 )a6(, n ) y (J3M 3 )(,nj - 2 - y Ml M 2 M 3 (J^Ol J 3 0) ' ^ 

rf 2 ^M:V fi ) y (^M 2 )a & (fi)^(.3A/3)(fi) = / ^ 2 ^SV ft ) y (^)a6(ft)^3M3)(fi) 

1 + (_!)*+*+* (Ji2J 2 ,-2|J 3 0) 

2 y Ml M 2 M 3 (j l0 j 2 0|J 3 0) ' 1 ' 

-Pf, V TB ab (f>\V TE (r\\Y (f\\ - i 1 ~ (-l) Jl + J2 + J3 r J 1 J 2 J 3 (Ji2J 2 , — 2| J3O) 
1 nY (J 1 AI 1 )( n ) Y (J 2 M 2 )ab{ n ) Y (J 3 M 3 )W - -I ^ yM 1 M 2 M 3 (J l0 J 2 0|J 3 Q) ' ( 66 - ) 

We now use Eq. (94) of Ref. [I?}, along with the orthonormality conditions, 

yL (AW VE/VB,ab ( ^_ (] yh ,~ w TE/TB,ab ( ~)_ () yVE/VB ,~.yyTE/TB,ab/*.\ « ffi ^ 
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to find that for J\ + J 2 + J 3 = even, 

,T,TB,ki („\ir,TB,k2,abf„\, Tr k 



Mt^ )ab (^ I r*drjj 3 (k 3 r) 



(Jl2J 2i -2|J 3 0) (TB,TB) (] s ATB,TB) ( , s 

(J!0J 2 0|J 3 0) Jjl >' l* ir ^,* ^ 



( JllJa, ~1| J3O) .{TB,VB) (f , .(TB,VB) , 

(JiOJ 2 o|j 3 o) ( * ir,J *. < (fc2r) 

(68) 



j3 .if TE ,k± 

ax w, . ' 



3 { lf L) ^r)jf 2 f L \k 2 r) 



(JllJ 2 , -1|J 3 0) {TE,VE) ( j , {TE,VE) , , , ( Jl2 J 2 , -2| J 3 0) (TE,TE) , , (TE,TE) } 

(W.O^O) ^ t* 1 ™*,* ^r) + (J l0 J 2 0|J 3 0) l* 1 ™*.* ^ 2r) 



(69) 



and for Ji + J 2 + J 3 = odd, 



»*™i)a6W*cffi) a0 ( x )*&,M,)W = ^ 3 m 3 / r^r^ 3 (fc 3 r) 



-.TEM.ab 



J1J2J3 



(JllJ 2 ,-l|J 3 0) .(TB,VB) x ATE,VE) ( , 

(JiOJ 2 0|J 3 0) ^ 



(Jl2J 2 ,-2|J 3 0) .(TB,TB) ( , s ATE,TE) x 

(Ji0J 2 0|J 3 0) ^ (*i r ^ a ,t 



(70) 



Here, we have used the radial functions defined in Eq. ([24"]) and Eq. ([125]) . Eq. ([58"]) . together with the overlap integrals 
in Eqs. ([55]). ([6U]). and ([701), provides the TB-TB-scalar, T^-TS-scalar, and TB-T E-sc&lai TAM-wave bispectra in 
terms of B t t s (ki, k 2 , k 3 ), the Fourier-space tensor-tensor-scalar bispectrum. 



VI. ONE SYMMETRIC TRACELESS TENSOR AND TWO SCALARS (OR TWO LONGITUDINAL 

VECTORS) 



We now calculate the three-point correlation that involves one symmetric traceless tensor and two scalars. The 
traceless tensor may be a longitudinal (L) mode, one of the two vector modes (VE or VB), or one of the two 
transverse-tensor (TE or TB) modes. The three-dimensional galaxy-survey observables of such a three-point tensor- 
scalar-scalar correlation were studied in Ref. [l5|. Moreover, a tensor-scalar-scalar correlation between inflationary 
gravitational waves (the tensor mode) and the primordial curvature perturbation (the scalar mode) arises generically 
during inflation fl3| . The TAM-wave bispectra we present here may be useful in calculating full-sky observables 
associated with these couplings. 

To be precise, we consider here a tensor-scalar-scalar bispectrum B sst (ki, k 2 , A; 3 ) defined by 

V^»(k 1 )VV»(k 2 )f cd (k 3 )) = (27r) 3 ^(k! +k 2 + k 3 )g ac g bd k 1 k 2 B sst (k 1 ,k 2 ,k 3 ). (71) 

Given that V a = V a </> is a longitudinal- vector field, we can also write this three-point function as a vector- vector-tensor 
bispectrum B vvt (k\, k 2 , k 3 ), where here "vector" is a longitudinal vector, of the form, 



y a (k 1 )y fe (k 2 )f cd (k 3 )) = (27T) 3 5 D (k 1 +k 2 + k 3 )g ac g bd B vvt (k 1 ,k 2 ,k 3 ). 
The three-point function of one TAM coefficient and two scalar TAM coefficients can be calculated from 



rL,k 2 



^a,k 3 



,ai3 



(72) 



(73) 



and using Eq. (fTTj) and A 



JM 



^ 1 m 1 )^^) t (Xm3))=-(4-) 3 h) 



That is, we only need to calculate the overlap, 



B sst (h, k 2 , fc 3 ) / ^^(xjtf^jtxjtf^k^fx). (74) 



L.k 2 b 



d A x^ 



(TOW* 



(J3M3 



,ab( X ) 



d A x 



(J2M2) 



(x) 



a,k 3 
(J 3 M S 



)„*(*). (75) 



for a = L, VE, VB, TE, TB. Below we discuss the three cases where the traceless tensor is longitudinal (L), vectorial 
(VE/VB), or transverse (TE/TB). 
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A. The transverse-traceless components 



We start with transverse-traceless tensors, a = TE, TB. Since the transverse-traceless fields are divergence-free, 
we can integrate by parts in the integral in Eq. (|75[) to get, 



1 



k\k 2 



We next use the definition [lj 



(76) 



(77) 



of longitudinal-tensor TAM waves, and the fact that the tensor mode is traceless, to rewrite the three-wavefunction 
overlap as, 



■ f L,k 2 b 



(78) 



for a — TE, TB. What we now need to calculate is the overlap of a transverse tensor, a longitudinal tensor, and a 
scalar. Using Eq. (|22|) and Eqs. (|62|) - (|66|) . we find that when Ji + J 2 + J3 = even, 



a X ^(J 1 M 1 ) 



x)*&X(x)*s: )Qb (x) 



1 



(j 3 M 3 )ab^J - ^9m?i£mAj^ I ^drjj, (fcir) jf 2 f } (k 2 r) jfJ' L) (for) 



{-hi^-h, — l|JlO) ,(i,VE) /, x -(TE,VE) 



J32J2, — 2| JjO) .(L,TE) s .(TE,TE) 



(J30J 2 0|J!0) (fc2r) ' 7j - t ( ^ ) + (J 3 0J 2 0|Jx0) 

and when Ji + J2 + J3 = odd, 



|(1h2) 



(79) 



(J 3 U 2 ,-l|JlO) .(£ 



Jt'^ Mjy 3 Z' v "'(for) + 



V6- 

(TB,VB) 



(J 3 2J 2 ,-2|JiO) ■ (L,TE) /, \ .(TB,TB 



(J 3 OJ 2 0|JlO) JJ2 ' t V /JJ3 '* V ' (J 3 0J 2 0|Jl0) 

Here, again, we use the radial functions defined in Eq. ([2^]) and Eq. (j2"5j) . 



j'L't (MA 



(M 



+ (lo2)|. 
(80) 



B. The vector components 

Next we consider the two vector modes, a — VE, VB. The TAM waves for the vector components of the symmetric 
trace-free tensor can be related to the TAM waves for transverse-vector fields 17[ through 



a*(./M) (X) = -^*(./A/)( X )> ^ a $ (JM) (X) = -_* (jM) ( X ). 



1 



(81) 



TE /TB k 

Since the (jM)ab' are traceless, we can do the overlap integral in Eq. (|75[) by integrating by parts to obtain 

'k 



^ 3 ^& 1 M 1 )a(x)*& 2 ) 6 (x)*&:;(x) = 



6 
1 

71 



,(82) 
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where a = TE, TB for tensor-valued TAM waves, and a = E,B for vector- valued TAM waves. Following the same 
techniques as in previous sections, we find for J\ + J 2 + J 3 =even, 



^-*f7^)a(-)*f/ 2 M 2 ) b ( X )*S3 3 r b ( X ) = V^^V 



,L.fc 2 



,VEM,abi 



r 2 drj Jl (kir) \ j£f ] (k 2 r) j^f ' L) (k 3 r) 



(J 3 lJ 2 ,-l|JlQ) (L,VE),, ,.(VE,VE) ( , s , (J 3 2J 2 ,~2|JiO) J L ,TE).(VE,TE) 



(J 3 OJ 2 0|JiO) 



(J 3 0J 2 0|.7i0) 



^g^A/ 3 M 2 Mi 



(J 3 lJ 2 ,-l|JlO) ,(L, 



+(1 2) 



(83) 



and for J x + J 2 + J 3 =odd, 



,L,k 2 



,VB,k 3! ab, 



(J3IJ2, — 1| JlO) .(L,VE) f , \ .(VB.VB 



■h?- -h, — 2| JjO) .(L,TE) ,{VB,TB) 



(JsO^OlJxO) (/c ™3, t + (j30j20 | Ji0) Jj 2 ,t Jj 3 ,t j + l 1 ^ 2 ) 



V8 



y M 3 M 2 Mi 



(84) 



Again, we have used the radial functions defined in Eq. (f2~i)) and Eq. (|2~5 



C. Longitudinal part 

Finally, we consider the longitudinal mode of the traceless tensor field. Once again, we can write the longitudinal 
tensor- valued TAM wave function ^ (jM) a b( x ) m t crms °f differential operators acting on scalar- valued TAM waves 
*(jm)( x )' as show 11 in Eq. (|77|). We can also put the two longitudinal vector- valued TAM wave functions into this 
operator form. Integrating by parts allows us to move around those gradient operators, and we re-write the overlap 
of three TAM waves as 

/ ^ 3 ^(^}:)(x)*Kii(x)*s /3)ab (x) =-L / {|*& lMl) (x)^-;(x)^ 3)afe ( X ) 

+ |*(} 1 M0( X )*&XW*SM3)a( X ) + a ^ 2 )} ■ ( 85 ) 

The first term is the overlap of two longitudinal tensors and one scalar, and the second term is the overlap of two 
longitudinal vectors and one scalar. The integrals are non-zero only if J\ + J 2 + J 3 = even, and the end result is 

/ ^oWCA^Kia^) = ^mXm, {/ r 2 dn. h (hr) [| (jf 2 f {k 2 r)jf 3 f (far) 

(J 3 lJ a ,-llJiO) .(l,ve) ( , x (l,ve) x (J 3 2J 2 ,-2|J 1 0) (L[T£) (iiTB) 
<J 3 OJ 2 0|JiO) ^ (fc2r) ^3, t ( fc 3?)+ (j30 j 20 | Jl0) J Al t (fer)j j3 , t (fc 3 r, 



I fa? Mj? 3 ? (far) - {J $£^ ti? Mtf 3 f (far; 



+ (1 o 2) , (86) 



_k 3 
k 

where the various radial functions have been defined in Eq. (|24| in Eq. (|25[) . 

We have thus calculated the overlaps of a traceless-tensor TAM wave and two scalar (or longitudinal-vector) TAM 
waves for the longitudinal, vector, and transverse-tensor components. These results, when inserted into Eq. (|73|) or 
Eq. (|73|) . give the TAM bispectra for a traceless tensor field with either two longitudinal fields or two scalar fields. 
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VII. CONCLUSIONS 

In this paper, we have calculated several bispectra for scalar, vector, and tensor fields in the total-angular- 
momentum formalism. We began with the scalar-scalar-scalar bispectrum. We then considered an example of a 
vector-vector-scalar (where here vector is a transverse vector) bispectrum and a tensor-tensor-scalar (where here 
tensor is a transverse-traceless tensor). The vector- vector-scalar bispectrum is of the form that may arise from the 
correlation of a magnetic field with a scalar field [ll|, while the tensor-tensor-scalar correlation is precisely the same 
form that arises in inflation (T^]. We obtained bispectra for TAM waves in the E/B basis (for vector fields) and the 
TE/TB basis (for the tensor fields) through intermediate steps that involved the TAM-wave helicity basis. We then 
moved on to calculate the TAM-wave three-point function for a tensor-scalar-scalar correlation that comes from a 
tensor-scalar-scalar bispectrum of precisely the same form that arises in inflation [l3| . For completeness (and to follow 
through in Ref. [25[ on CMB signatures of correlations of the form considered in Ref. [HI]): we have also considered 
the bispectra for two scalars and either the longitudinal of vector components of the tracclcss tensor field. 



type of correlation 


Fourier-space bispectrum 


TAM-wave result 


scalar-scalar-scalar 


m 




vector-vector-scalar (transverse vectors) 


m 


(S3 and dSD-dMJ 


tensor(T)-tensor(T)-scalar 


(23 


(58) and (|68)-(70) 


tensor(T)-scalar-scalar 


CD 


(73). (79), and (80) 


tensor(T)-vector-vector (longitudinal vectors) 


(72) 


(73). (79), and (80) 


tensor (V)-scalar-scalar 


mo 


(74). (83), and (84) 


tensor(V)-vector-vector (longitudinal vectors) 


(72) 


(73), (83), and (84) 


tensor(L)-scalar-scalar 


CD 


(74) and (86) 


tensor(L)-vector-vector (longitudinal vectors) 


ca 


(73) and (86) 



TABLE I: List of the types of three-point functions we consider, the equation where the Fourier-space bispectrum we consider 
is defined, and the equations that contain the central results for the TAM-wave three-point functions. We distinguish between 
longitudinal and transverse vectors. The labels T, V and L refer to transverse tensorial, transverse vectorial and longitudinal 
part of a traceless tensor field respectively. 

In the plane-wave formalism, the three-point correlation in Fourier space is parametrized by a bispectrum 
B(ki,k2,k 3 ) that depends on the three wavenumber magnitudes ki,k 2 ,k 3 only, multiplied by a "momentum- 
conserving" Dirac delta function, a consequence of statistical homogeneity. In the TAM formalism, generically the 
three-point correlation is parametrized by exactly the same bispectrum, but multiplied by an "angular-momentum- 
conserving" Clebsch-Gordan coefficient (which we write in shorthand as a Gaunt factor), a consequence of statistical 
isotropy (the Wigner-Eckart theorem), and an integral over radial profiles which does not depend on the azimuthal 
quantum numbers. Statistical homogeneity along the radial direction is encoded in the integral over the radial profiles. 

Since the E/B and TE/TB TAM waves have definite parity, parity conservation is manifest in the bispectrum 
(determined by whether J\ + Ji + J3 is even or odd). This is to be contrasted with the plane- wave formalism, where 
the two linear polarizations of a transverse-vector or a transverse-traceless-tensor are identified with each other by 
a simple rotation about the direction of wave propagation. In the TAM formalism, the E/B modes have different 
geometrical factors and different radial integrals in the bispectrum, despite the fact that their power spectra must 
coincide by statistical isotropy [l?}- This means that E/B modes correlate to other fields differently when cubic 
interactions are taken into consideration. They thus have different bispectra. 

Our survey of bispectra was not at all exhaustive. We showed in Section ITVAl for example, that vector- vector-scalar 
correlations can be parametrized, assuming parity invariance, in terms of two bispectra, and we then considered one of 
those. There may likewise be additional types of bispectra involving tensor fields and beyond those that we considered 
here. There are also additional forms for bispectra that may arise if we allow for parity violation, and in this case, 
the helicity-basis TAM waves we developed in Ref. should be particularly appropriate. 

The value of the TAM-wave bispectra we have discussed here becomes apparent when we realize that many cos- 
mological measurements are performed on the full sky, or over wide-angle surveys. Angular correlations are then 
decomposed into angular power spectra Cj (usually written as C/) parametrized by multipoles J. The advantage of 
TAM waves is that once the proper scalar, vector, or tensor TAM waves are identified, rotational symmetry guarantees 
that the observable Cj will receive contributions only from TAM waves with that same J. In this sense, TAM waves 
provide a more natural choice of basis functions than the conventionally used plane waves. The advantage for three- 
point functions is that once the proper TAM waves are identified, the results we have derived in this paper provide the 
TAM-wave angular bispectra in terms of the more commonly seen Fourier-space bispectra. Since observable angular 
bispectra will be obtained from projections of the TAM-wave bispectra, the Ji, Mi, J2, M 2 , J3, M 3 dependences of the 
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angular bispcctra can be read off directly from the results we have presented here. For calculating CMB observables 
from bispectra involving vector or tensor fields, the TAM approach is particularly powerful in that equivalent results 
can be obtained without any Wigncr-6j or Wigner-9j symbol, nor does any additional summation over the orbital 
angular momentum arise, as opposed to using plane waves. 

One example of the utility of this formalism will be provided in forthcoming work (25j where we calculate the 
bipolar power spectrum of the cosmic microwave background that arises from scalar, vector, and tensor distortions 



to the local density two-point autocorrelation function of the form discussed in Rcf. 15]. In other work [26f we will 
show how these TAM-wavc bispectra can be used to construct estimators from wide-angle galaxy surveys (including 
rcdshift-space distortions) for the types of distortions considered in Ref. fl5j . 
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Appendix A: TAM-wave overlap integrals in the helicity basis 

In this Section, we calculate the overlap of three TAM waves in the helicity basis. We first calculate the overlap of 
two transverse- vector TAM waves and one scalar TAM wave in Section IA II and then the overlap of two transverse- 
tensor TAM waves and one scalar TAM wave in Section IA 21 

Our starting point will be the overlap [27j . 

rl 2 A v, (A\ v, , (ft\ v, (A\ — ( y\h+h+h+ nhhh {hs\hs2\h - S3) 
a n sl y( iimi )(nj S2 r(; 2m2 )(nJ S gY( i3TO 3)(nJ — [-L) y mi m 2 m 3 — JhQj^m^) — ' 

of three spin-s spherical harmonics, which holds for s\ + S2 + S3 = 0. 

1. Two transverse vector TAM- waves and one scalar TAM wave 

Let us first calculate 

' d 2 «-9 ab Y ( X j[ Ml)a (n) Y{/ 2M2)b (A) Y (JsM3) (A) , (A2) 
for Ai, A2 = 0, ±1. We use the completeness relation for the spherical basis, 

g ab = 5>l(A) £ f (A) = £ (-l) A ^(n)e fc _ A (n). (A3) 

A A 

The g ab on the far left does not really depend on n, but the decomposition can be done at each h. Then, we find 

/ rf 2 ft.9 ah ^ 1 1 M 1 )a(ft)^M 2)6 (ft)^(./3M3)(ft) = E(-!) A / rf2 ^A(fi)^A(fi)^ ( A / 1 M l)Q (ft)^M 2 ) fc (ft)^(^M3)(ft) 

= (-l) Al S Xl -x 2 J d 2 h- Xl Y { j lMl) {h) Xl Y { j 2M2} (n)Y { j 3M3) (h). (A4) 

Here, we use Eq. (52) of Ref. [l7| to relate the vector spherical harmonics in the helicity basis to the spin-weighted 
spherical harmonics [18| 

^,(n)Y( A JM)a (n) = _ A Y (JM) (n)<W. (A5) 

Since the sum of the three spins in the three spin-weighted spherical harmonics in Eq. (|A4[) is zero, we use Eq. (|A1[) 
to calculate the angular integral of three helicity-basis TAM waves to be, 



/ 



^% a V 1 M 1 )a( fi ) r ( A /:M 2)b (fi)^(^M 3 )(n) = S Xu . X ^Xm 3 <Jl ( j'oj 20 | jffi 0) ■ ( A6 ) 
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Due to the Kroncckcr delta, the only non-zero combinations are 

|d 2 % afc r ( ° /iMi)a (n)y ( 72M2)b (n)r ( , /3M3) (n) = g££ M „ 

J d n 9 r (J lMl )aWy(J 2 M 2 )bWHJ 3 M 3 )W - -y MlM2M3 (j^j^j^ ' 

/ d 2 n 9 ^ Mi) Jfi)Y ( y 2M2)b (A) Y (JsM3) (h) = ^q)^^ ' (A?) 

2. Two transverse tensor TAM-waves and one scalar TAM wave 

We now calculate the helicity-basis overlap integral, 

-) 2 A Tab,cd v Xi f&\V^2 



^^ O0,Ca ^U 1 )a6(ft)^M 2 )c d (fi)^(7 3 M 3 )(fi), (A8) 

for Ai, A2 = 0, ±1, ±2. With the completeness relation, 

iab , cd = STif(h)if*(n) = (-l) A ef (fi)e-A(fi), (A9) 

A A 

and Eq. (100) of Ref. [13] (the relation between tensor spherical harmonics and the spin-weight s = 2 spherical 
harmonics) , 

i<$(n)Y ( x JM)ab (n) = _ A F (JM) (n)5 AV , (A10) 

the integration becomes 



A J 

= (-!) Al <*Ai,-A a / rf 2 n_ Al Y ( j 1 A/ 1 )( fi )A i y(J2A/2)(n)y(j3M3)(n). 



Then we find 



(All) 

d 2 n^/;- ) (n)r ( ^ 2M2)ab (n)r, /3M3) (n) = (-1)^,-^^ ^^'q^ • ( Al2 ) 



Appendix B: A semi-classical interpretation 



In this Appendix, we briefly discuss a semi-classical picture which allows us to interpret our results for the three- 
point function involving two transverse vectors and one scalar in the limit of large angular momentum. 

In Sec. IIV1 we presented three-point correlations involving two transverse vectors and one scalar, namely Eq. (|47p . 
together with the overlap integrals in Eqs. (f5T|) . (|52"1) . and ([5U|) . Some re-arrangment of those results enables us to 
rewrite them as 



( v (fX)\fXAw) = ~ ( 4 -) 3 H) Jl+j2+j3 b vvs k 2 , k 3 ) g( }: 



M 2 M 3 



x / r dr 



[1] + [2] - [3] , EtE) {EiE) „/mToT J ' Jl (fcir) JJa (fc2r) 



(Bl) 



( t/ (j;M 1 ) y (./2A/ 2 )^(.73M 3 )) = ( 47 (-*) S ™ fc 2, fe 3 ) 2 y[Xp g ^i^3 

x / r 2 drj. h (k ir ) jj 2 (k 2 r) j, h (k 3 r) , (B2) 
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\ ViMO V,M a ) VsMi)' _ W ^' fc2 ' ^ M1M2M3 (J l0 J 2 0|J 3 0) V 2J 3 -1 

< V(- Ji + J 2 + j 3 ) (Ji -J2 + J 3 ) (Jl + J 2 ~ J 3 + 1) (Jl + J 2 + J 3 + 1) 

x / r 2 d nJl (for) 3 fJ ] (k 2 r) jj 3 (for) . (B3) 



Here we have introduced the shorthand notation [i] = Jj(Jj + 1) for i = 1, 2, 3. These three-point functions have now 
been put into forms comparable to the three-point function, Eq. (|38p. for three scalars. We have a bispectrum function 
in terms of wave vectors with some suitable multiplying prefactor, a Gaunt integral in accord with the Wigncr-Eckart 
theorem, and an integral of three radial profiles. For a vector-vector-scalar correlation, we have an extra geometric 
factor made of angular momenta Ji, J2, and J3. In the limit of large angular momenta, i.e. J%, J2, J3 S> 1, these 
factors are reduced to cosine and sine 



[1] + [2] - [3] y/(- Jl + J 2 + J 3 ) (Jl - J 2 + J 3 ) (Jl +J2-J3 + 1) (Jl +J2 + J3 + 1) . a m „. 

, — ~ cos(7i9, , ~ smtfio. (o4) 

2VTTTI2T 2yiipl 

where 812 is the angle between J\ and J2, in a triangle whose three sides have length Ji^Ji and J3; rcpectivcly. 

These geometric factors, interpreted in a semi-classical way, can be attributed to the vector structure of vector 
TAM waves. First consider the case of two B-mode transverse vectors. The three-point function is proportional to 
the overlap, 

rf 3 ^(};f;:)(x)*^ 2)a (x)*^ M3) (x) = / d * x (^*!W*)) • (^*!W X >) *<W*) 

1 f * 



' T '(J 3 M 3 ) 



= VW\V¥}c°s0i2 J d 3 xy { j lMl) y { j 2M2) y 



= -COS#i2 J d 3 X^ {JlMl) ^ { j 2 M 2 )^(J 3 M 3 ), (B5) 

where according to the construction in Ref. [l7| we have used ^^-.(x) = K a ^ JM ^ (x) /y/ J(J + 1), with K a = 

—iL a and L a being the orbital angular-momentum operator. In the second line, we associate two classical angular- 
momentum vectors Ji, J2, of magnitude J\(Ji + 1) and J2(</2 + 1) respectively, with the first two TAM waves. Note 
that Ji and J2 are not fixed vectors, since they precess about the z axis. But by angular-momentum conservation 
they always differ by a third angular-momentum vector J3 of magnitude Js(Js + 1) that is associated with the third 
TAM wave, which also precesses about the z axis, as shown in Fig. [TJ With this picture, the cosine factor then arises 
naturally. 

The case of one B mode and one E mode can be analyzed similarly, 



^^(};Z)(x)*(};m 2)q (x)*^ M3)( x) = / d* x 1 — *« }iMi)( x) 1 • 1 ^(x) 1 * (j3A/3)( x) 



1 



d 3 .TJi* (JlA/l) • (k 2 x J 2 *(./ 2 a/ 2 ) ) * 



= (Ji • fk 2 XJ 2 )1 / d 3 x^ (JlMl )^(j 2 M 2 )^ 



(J1M1) ' X J2*(J 2 M 2 ) ^ *(J 3 M 3 ) 

(,/iM 1 )*(J 2 M 2 )*(J 3 Af 3 ) 



((ji x j 2 ) -k a ) /" d 3 x* (JlMl) * (j2M2) * (J 3 M3) , (B6) 



where this time we have also used ^(j M ) a ( x ) = Mj*^ (x) j %J J( J + 1 ) , with M a = ie a f )c V b K c /k. We clearly see a 
factor |Ji x Ji I = sin 0i2, which is the sine we find in the large- J limit. 
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FIG. 1: Geometry of the semi-classical picture: three classical angular-momentum vectors Ji, J2, and J3 form a triangle and 
precess about the chosen z-axis, associated with the three TAM waves of wave numbers k\ , k2, and £3 respectively. The angle 
9i2 between Ji and J2, however, remains the same. The wave vectors ki and k2, which are not definitive for TAM waves, are 
perpendicular to angular momenta Ji and J2, respectively. 



To close, we consider the case of two E-mode vectors 

/ HSdW^)»!U)W = / d3x (^*&^oW 



V x L . 



*(} 2 M 2 )( X ) *(J3M 3 )(X) 



1 



d A X (ki X Ji ) *( JlMl) • (k 2 X J 2 *(J 2 jU 2 ) J *(J 3 J\/ 3 ) 

#1 



ki x JiJ • \ k 2 x J 2 J / d a;*(j 1 Mi)*(j 2 M2)*(J3M3) 

ki ■ k 2 ) (Ji ■ J 2 ) - (ki • J 2 ) (Ji ■ k 2 ) J d 3 x^ { j lMl) ^ { j 2M2} ^ { , hM3} . 

(B7) 

We then recognize a term proportional to Ji • J 2 , i.e. cos# 12 , and a second term which has no such factor. It closely 
resembles the structure in the exact result Eq. (|B1[) . 

So far, we have restricted our discussion to the three-point function involving two transverse vectors and one scalar. 
Besides, the semi-classical picture we have proposed is still insufficient to quantitatively pin down the correct forms of 
the radial integral and the Ji, J 2 , J3 dependence. However, we may gain insight from this picture for other three-point 
functions, such as the ones involving transverse-tensor field. 
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